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RESEARCH MEMORANDUM 


A METHOD FOR THE DKIIGN OF SWEPTBACK WINGS WARPED TO 
PRODUCE SPECIFIED FLIGHT CHARACTERISTICS 
AT SUPERSONIC SPEEDS 
By Warren A. Tucker 

SUMMARY 


One of the problems connected with the swepthack wing is the diffi- 
culty of controlling the location of the center of pressure and hence 
the pitching moment. A method is presented for designing a wing to be 
self -trimming at a given set of flight conditions. Concurrently, the 
Bpanwise distribution of load on the wing is made to be approximately 
elliptical, in an effort to maintain low wing drag. 

These flight characteristics are achieved by warping the wing out 
of a plane. The required warp is determined by the values of the coef- 
ficients of a four-term series describing the pressure distribution; 
these values in turn are determined from four conditions on the lift, 
pitching moment, and spanwise load distribution. 

The method is dir ectly applicable to several wing plan forms, 
including the triangle and the sweptback plan form with finite tips, 
under the restriction that the leading edge must be subsonic and the 
trailing edge supersonic. The application to any specific problem is 
simplified to a routine computational procedure by the presentation of 
certain basic data in tabular form. A discussion is given of seme 
points to be considered in the application of the method to a practical 
case, and several representative examples are worked out. The resulting 
wings are shown to be ones which might practicably be built. 


mERODUUTIGN 


The evolution of the sweptback wing for efficient flight at super- 
sonic speeds has reached the point where the stability and control prob- 
lems are being investigated. This situation implies that not only the 
lift and drag of the wing but also the pitching moment must be considered 


in relation to the airplane as a whole. 
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In order to be truly efficient at the design Mach number, the wing 
should produce the design lift coefficient without creating about the 
airplane center of gravity a pitching moment. that would require a large 
deflection of the trimming device (with a correspondingly large drag). 

In addition, it is generally desirable that the spanwise distribution of 
lift be as nearly elliptical as possible and that any adverse pressure- 
gradients on. the wing be small so as to retard separation of the flow. 
These two conditions are not sufficient to guarantee that the wing drag 
will be a minimum, because at supersonic speeds the drag due to lift is 
also dependent on the chordvise loading; they are, however, conducive 
to low wing drag. 

The use of wings warped to produce a constant pressure over the 
surface has been proposed to eliminate the large adverse pressure gradi- 
ents encountered with the flat wing. For a given plan form, however, a 
uniform pressure distribution allows no control over the pitching moment. 
The wing warp necessary to produce certain other pressure distributions 
has been derived (reference 1), but these distributions do not lend them- 
selves readily to the control of pitching moment; in fact, the conical 
nature of the pressure distributions fixes the center of pressure at the 
center of area for triangular wings. 

In the present paper, data are presented from which the wing warp 
necessary to produce a certain type of pressure distribution may be 
determined. A development is then given In which certain constants 
appearing in the expression for the pressure distribution are determined 
by conditions on the lift, pitching moment, and spanwise load distri- 
bution. In this manner a method is derived for designing a wing of given 
plan form, operating at a given supersonic Mach number,__ to have a speci- 
fied lift coefficient, a specified center of pressure, and a nearly 
elliptical spanwise load distribution. Although the pressure gradients 
are not controlled directly in the method, the type of pressure distri- 
bution used ensures that for most reasonable design conditions the 
gradients will not be excessive. , There is no reason to believe that a 
configuration using a self-trimming device designed by this method will 
necessarily have a lower drag than will a similar configuration using a 
flat wing and a deflected trimming device. The possibility does exist, 
however, and should probably be investigated. ; 

The method is applicable to a wide class of wing plan forms shown 
in figure 1; the principal requirement Is that the leading edge must be 
subsonic and the trailing edge supersonic. The presentation is made in 
a form suitable for engineering use, and a table and computational form 
are provided so that the application of the method Is reduced to routine 
computation. — ■ 
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SYMBOIB 


c( 1 - \) ■ 



■"O = 

(1 

- k) [l - ff (l - \ )] 



- 1 


C 

local chord 


c r 

root chord 


c 

mean aerodynamic chord \— f* c 2 dy = ^{?- + ). <•. \ 

V s Jo 3d + x) 


C D 

drag coefficient 


C Z 

lnrfl.1 lm. rw»-pfir.-f«n+. (Lift on chordwiae strip dy in width) 

* 

V qc / 


C L 

lift coefficient 

l IS ) 


Pm 

pitching-moment coefficient, positive when pitching moment 



^ , ^Pitching moments 

tends to move wing apex up ( 1 

\ qSc / 


ksl 

“l 




taper ratio c ^ ord - 

\Root chord/ 


m 

cotangent of sweephack angle of leading edge (see fig. 2) 


m l 

cotangent of sweephack angle of trailing edge (see fig. 2) 


M 

free- stream Mach number 


ns pm ( 

see fig. 3) 


P L 

static pressure on lower surface of wing 

■rf 

% 

static pressure on upper surface of wing 
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P lifting-pressure coefficient 

q free- stream dynamic pressure 

r e (see fig. 3) 
s semispan (see fig. 2) 





S 


wing area 


x,y rectangular coordinates parallel 'and normal, respectively, to 

free stream, with origin at wing apex (Bee fig. 2) 

x ' distance behind leading edge; measured in free- stream 

direction - : - ■ ~ 

Xq distance of moment axis behind wing apex (see fig. 5 ) 

x 7 distance of moment axis behind leading edge of mean aero- 
dynamic chord (see fig. 5) ! 

z distance perpendicular to xy-plane, positive up 


ANALYSIS 

General 


A convenient method is derived in reference 2 for finding the wing 
shape corresponding to a given pressure distribution. In the present 
paper, the lifting-pressure distribution over the wing is taken to be 
of the form I * 


P = Ci 1 + C 2 'x + C 3 '|y| + C^'y 2 (l) 

where the axis system is that shown in figure* 2, and C 1 *, C 2 ’, C^* , 

and C^' . are as yet arbitrary constants. Other terms could have been 

included in the series but the terms shown gave acceptable results with- 
out requiring undue labor. For convenience, the coefficients of the 
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series may be replaced by others similar in nature Buch that the lifting- 
pressure distribution can be expressed by the following equation: 


_P_ 

C L 


SL + 1 - ^ £2 _x_ . c 3 

C L 1 - X °L c r C L 



( 2 ) 


For purposes of calculation, the wing is assumed to have no thicks 
ness so that the shape derived is actually the mean surface of the wing. 
Within the assumptions of the linearized theory, an arbitrary thickness 
distribution, symmetrical above and below the mean surface, can then be 
added with no effect on the lift and pitching moment. 

Suitable integrations of the pressure distribution over the plan 
form may be performed to obtain equations for the lift coefficient, the 
pitching-moment coefficient, and the spanwise load distribution. One 
condition may then be imposed on the lift, one on the pitching moment, 
and two on the spanwise load distribution. This procedure results in 
four linear equations in the four unknowns C^J C L , C 2 /C L , c 3/ c l, 

and The values for these constants may then be substituted into 

equation (2), and the shape of the wing (that is, the warp) necessary to 
produce this pressure distribution can be found by the method of 
reference 2. 


The foregoing material has described the method in general terms. 

In the following sections more detailed descriptions are given: First, 

of the procedure used to find the warp corresponding to each component 
of the pressure distribution; second, of the method used to determine 
the constants for the case of plan forms having pointed tips; and last, 
of the corresponding procedure for plan forms having finite tip chords. 
Although the determination of the constants is in principle the same for 
both types of plan forms, certain simplifications occurring for the 
pointed-tip case make not only the actual numerical work of determining, 
the constants but also the exposition of the procedure simpler for this 
case. 


Warp 

The warps or wing Bhapes necessary to produce the several components 
of the pressure distribution given in equation (2) are first found sepa- 
rately as functions of" the four constants Cg/C L , 

and Later, after numerical values of the constants have been 

determined by the conditions Imposed on the lift, pitching moment, and 
spanwise load distribution, the separate shapes Eire superimposed to form 
the finEl war^ped wing. 


1 imuuui iUP^' 
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The general idea In finding each wing shape Is first to determine 
the slope of the wing surface q associated with the pressure 

distribution under consideration and then to integrate this slope in 
the x-direction to get the z-ordinate at any' point (the direction of z 
is taken mutually perpendicular to x and y, positive upwards). Of 
the available methods for finding the wing slope corresponding to a 
given pressure distribution, that presented in reference 2 was chOBen 
for the particular problem. The principal advantage of this method is 
that it eliminates the need for considering z in the integrations 
involved and so simplifies the integrations. ; 

The slope of the wing surface corresponding to each term of equa- 
tion (2) is found by application of equations (8) and (17) of refer- 
ence 2 . The wing shape as a displacement from the z =_ 0 plane 1 b then 
found by integrating the slope in the x-direction; thuB, 



The following equations result for the wing shapes corresponding to the 
four components of the pressure: 


_° l © 


m 


xR 1 = 


m 


4 n 


2 'Jl - n 2 !- 2 - 2 cosh" 1 


_1_ 

nr 


Vl - n 2 (l + r)cosh _1 


1 + n 2 r 
n(l + r) 


+ 


Vl - n 2 (l - rjcosh"?- 


1 - n 2 r ~ 
n(l - r) 


( 4 a) 




x 2 ^ \Jl - n 2 !* 2 
■ 4rt 


- 2r 2 cosh“^ - 



+ 



n2 d 


cosh" 


1 + n 
n(l + r) 


+ 


1 

■ *» 
n 2 (l - r 2 ) 2 

cosh"^- 

1 - n 2 * 


l/l - n 2 

2 r 

n(l - r) 




( 4 b) 
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The significance of the quantities r and n Is most clearly 
seen "by reference to figure 3. Calculations have been made of the 
quantities R^, Rg, and R^, which are in a certain sense the 


conical parts of the wing shapes, and the results are presented in fig- 
ure 4 and. table I. The figure is intended to be merely illu strative; 




8 


HACA RM L51F08 


study of the figure provides a qualitative Idea of the various wing shapes 
One Interesting fact to notice is that no infinities occur at the center 
line (r ■ 0) for the cases in which the pressure is proportional to x 
and to |yj (compare with the shapes derived in reference l) . 


Evaluation of Constants for Pointed-Tip Wings (X » 0) 

The pressure distribution given by equation (2) is integrated over 
the plan form to obtain. an expression for the lift of the wing. If the 
limits shown in figure 5 are used, the following equation expresses the 
value of the lift: 



If the 'indicated operations are carried out and the lift coef- 
ficient is formed, the following equation results: 


2 - t £2 . 1 £2 + 1 2fc 
3 C L * 3 + Z 


(6) 


The pitching moment about an axis a distance x Q behind the wing 
apex may also be found by the following equation: 


^ = -20, 




(x - x 0 )JL dx 
"L 


(7) 


After the pitching-moment coefficient is formed, the following equation 
is obtained: 


Pm 

C L 



1 Xq 

2 cr 


2 - k ~ 1 £l 

2(1 - k)J C L 
3 - k ~ [ c 3 

' 8(1 - k)J C L 



3 - 3k + k 2 
4(1 - k) 



14 - k £4 

20(1 - k) C L 


£2 

Cl 


+ 


( 8 ) 
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The spanwise load distribution is found from the following Integration: 


cc t 

°r c L 


if 

Cr Jy/m ^ 


where c is the local chord and Cj is the local lift coefficient. 
The integration has already been made in finding the lift. The 
following equation results: 


ffL.. (Si + k=-£.£a\ _ 

c r c L V°L 2 C V 




1 + k c 2 ^3 
V 2 ” C L °L 



o 2 -^ 

C L 


( 9 ) 


For purposes of reference, the spanwise load distribution if elliptical 
would be given by the following equation: 


ih- - SGLLiii/rr? do) 

c r C L * 

where 1=0 for the pointed-tip case. Equations (6), (8), and (9) are 
now used to find values for 0^0^, Cg /C l* an ^ ^kl^L' ® ie 

following conditions are first applied to equation ( 9 ): 



2(1 + 1 ) 
n 


(value for ellipse) 



(value for ellipse) 


r f 11 ) 


where again 1 = 0 for the pointed-tip case. These conditions are not 
quite arbitrary but were chosen after trial of a number of possibilities. 
The selection of these particular conditions not only made possible a 
solution for the four unknowns but also resulted in a single equation 
for the spanwise load distribution which was a fair approximation to an 
ellipse. (The, degree of the approximation is shown subsequently. ) 
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Substitution of 
values for 


£g m 4 2 £2 

Cl (l + k)rt 1 + k 


The solution for the constants may now be made, 
equations (11) Into equation (9) gives the following 
and . j 

£l = 4k ' 1 - k £2 

Cl " (1 + k)rt 1 + k C L 


These values are substituted into equation (6) and the following solu- 
tion for. Cj^Cl is got : 

c L * 


The values of 0-^0^ c 2f C L * C k/ C L 611,6 Bul3B ' bi ' fcu * e<i ircto equa- 
tion (8) and the left-hand side is set equal to zero to arrive at the 
solution for Cg^C L . The solutions are collected in the following 

equations : 


or 




• “ * 

°L * 


■ *r 

c 3 8(1 + k) 

3(4 - k) 

1 + 7 k + k 2 

Cl 1 - k 

10(1 - k) 

5k (l -■ k 2 ) 

C 3 8(1 + k) 
C L = 1 - k 

7 - 3k 
10 (l - k) 

1 + 7 k + k 2 
5« (1 - k 2 ) 

C 2 _ 4 

. ? °3 


°L (1 + k)« 

1 + k c L 

-= ■; 

C 1 4k 

+ 1 - k c 3 


C L. (1 + k)n 

1 + k C L 





> ( 12 ) 


The two forms for are given because in some cases the center of 

pressure is located more conveniently with respect to the' mean aerodynsmi 
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chord c. The geometrical relation between Xq and x' iB indicated 
in figure 5 ; the analytical relation is 


= 2 . il + 1 

c r 3 c 3(1 - k) 

For k = 0 (triangular wing) , equations (12) simplify to the 
following equations : ^ 

fit = 6 - * 


(13) 


£2 

cl 


\5 5* 2 c T ) 


or 


2l . e (J_ _ i . S') 

C L °\l0 5it 5 ) 


fa. it . 2 

-1 — 

J L 


> (1A) 


Equations (12) when substituted into equation ( 9 ) result in the 
following equation for the spanwise load distribution: 


cc z 

c rCL 


-* + C 6 -^) 02 - ( 6 -‘^)'’ 3 


(15) 


This equation is unif ormly valid for al i pointed-tip wings , independent 
of the lift, center of pressure, Mach number, and relative sweeps of the 
leading and trailing edges. The load distribution given by equation (15) 
is compared with the elliptical distribution of equation (10) in fig- 
ure 6. As a matter of incidental interest, the spauwise center of load 
on one wing panel is located 0.409 semispan outboard of the wing center 
line for the load represented by equation ( 15 ) , compared with a corre- 
sponding value of 0.424 for the elliptical load. 

Evaluation of Constants for Wings with Finite Taper 

The problem of wings with finite taper can be approached in two 
wayB. The more obvious method is to assume that the pressure distribution 
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defined "by equation (2) applies over the entire wing surface and to 
calculate the required warp. (A separate calculation for the warp in the 
tip region, which is shown shaded in fig. 2, is required, hut this cal- . 
culation is not impossible to make.) The disadvantage of this procedure 
is that at the tip the required wing slope takes on very large values 
(theoretically infinite). A more practical approach, and the one 
adopted in this paper, is to relax the condition on the pressure in the 
tip region. For a flat lifting wing with subsonic leading edges, the 
average pressure in the tip region is known to be close to zero (ref- 
erences 3 and ^) . It is not entirely illogical, to suppose that for a 
slightly warped wing the pressure in the tip will also be very small. 

If for the warped wing the pressure in the tip were taken to be exactly 
zero, the equations for the lift, pitching moment, and so forth would 
be derived by first integrating the pres sure. _distribut ion defined by 
equation (2) over the entire wing, including the tip region, after which 
the Integral of the same pressure over the tip region would be sub- 
tracted. In order to keep the equations within reasonable limits, a 
constant pressure was Instead subtracted from the tip region. The value 
of this pressure was taken to be the value given by equation (2) at the 
middle of the tip chord; namely. 


P_ 

Cl 




C 2 £3 . C 4 

°L C L C L 


(16) 


After uhe foregoing assumption has been made, the determination of 
the four constants proceeds very similarly to that for the pointed-tip 
wing: The pressure distribution is integrated to obtain equations for 

the lift, pitching moment, and spanwise load distribution; two conditions 
are set' on the spanwise load distribution, and the resulting two equa- 
tions are solved together with the equations for the lift and the 
pitching moment to give values for the constants C i/C L , CqJ C l , C^C l , 

and C^C L . The work Just outlined is now given. 


The lift equation, corresponding to equation (6) for the pointed- 
tip wing, is 


1 - (1 - A) 


t . 2 - k l g (l + k) 2 - 1(1 + k) . c 2 J 
[3(1 - 1) " 3(1 . x 2 ) ' 2(1 - 1) A J C L + 


1 + 21 , 


1 + 3* ■ A 

C|f 

3(1 + 1) A 

c L + 

6(1+1) A 

c L 


(IT) 
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where 

A= X 2 (l ~ k) 

(k + n)(l - X2) 


(18) 


In this equation, all the terms containing A arise from Integration 
of the pressure over the tip region, whereas the remainder of the terms 
represent the Integration of equation (2) over the whole wing area. If 
X Is set equal to Zero, the equation reduces to the form given in equa- 
tion (6). 


The pitching-moment equation, corresponding to equation (8) for 
the polnted-tlp wing, is 


(1 + X + X 2 


, c m [3(1 + X) >0 [3(1 - kX) (1 + k)(l - X) J 

1 2 “ B * 2(1 I k) " 2T1 - k) 

L J L 


-C 


3(1 - kX) (1 + k) (l - X) 
2(1 - X) ‘ 2 


X(1 - k) ♦ _2(1 - X) 

2(1 - xl 


B 


3sl Jl-JE + & + 


11 - X 


1^(1 - X) _ (1 + k + k?)(l - X) 2 
1 - k Ml - k) 



4 
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where 




« = 3X 2 (l - k) 

3 2(k + n)(l - XT 


c - 31 2 , l 3 (2n + k)(l - k)_ 

2(k + n) 2 (k n) 2 (l - l) 


> 


( 20 ) 


The remarks following equation (17), with A replaced by B and 
C, and equation (6) replaced by equation (8), also apply to equation (19) 

Because of the existence of a tip region for the wing of finite 
taper, the spanwise load distribution must be described by two equations, 
one applicable from the center line out to the beginning of the tip 
region and the other applicable over the r em ainder of the span: 


For 0 =1? a =? 1 


X(1 - k) 

(k + n)(l 1 \)’ 


cc i /fcji i - k 

<h?L = \Cl 2(1 - l) C L J 



(1 4- k)(l - 1) 

2 C L 


(1 




(1 - X) ^ cr 3 (21a) 
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F°r i - . %iL £o £ 1, 

(k + n)(l - X,) 


cc l 

c r°L 


J (1 + n)(l - X) C 1 , 

' i - k 


(1 - k) 2 - fe - 1(1 + k)11~l(l + n) - (k + nl] £2 
2(1- k) (1- A) C l 


X(l + n) - (k + h) 

• 


(1 + n)(l - X) C 1 

1 - k 

\CL c £/ 

* “ * 

1 - k c L 

r- 

— 




(1 - A) fek(l - k) - 2(k + n) + 1(1 + k)(k + n)1 £2 

2(1 - k)(l - A) Cl " 


(1 - k) - (1 - X)(k + n) (l - l)(k + n) £4l _ 

1 - k C L 1 - * C L f 


(1 + k)(i - x)_ °2 ( _ x) £i3 

2 c L + u c l 


,2 - (1 - 1) ^ ff 3 


£4 

3 l 


(21b) 


The conditions given in equation (U) are now applied to equa- 
tion (21a) to obtain the following relatione: 


£l 

°L 


4k(l + X) 
(1 + kT 5T 


+ 


1 - k c 3 
(1 - X)(l + k) ^ 


( 22 ) 


£2 4(1 - l g ) _ 2 2l 

Cl * (1 + k)rt 1 + k Cj^ 


( 23 ) 


These values of C^/Cl 
tion (17) and equation (19). 


and Cg/CL are substituted into equa- 
The center of pressure is fixed at Xq 
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Isy putting C m equal to zero in equation (19) and the following 
equations for and sure obtained: 


£4 _ 6(1 + X) 

C L fa. + 3X) - 6(l + X)A 


1 


8(1 + X) - 4x 2 
3« 


(4 - 2X) (l + X) 

it 


A 


(24) 


(1 - k)(l + X + X 2 + x3) c 3 J 
8(1 - XTTTTkJ L l£ 




11 4. 3fli - > 4 ■ 3 .( 1 + *1 4 . * ). g . e 

1 2 . 4(1 - k) 10(1 - k) b 


.£4 

°L 


. a(l.±X).J 

(l + k)it ' 


3! 


3(1 + k) .- (1 + k)(l - X) 2 - (2 - X)(l + k)B 
(3 + k - k g )(l + X) - (3 + 5k + k g )X g + (l + k + k g )x3 


*0 


■J 


2(1 - k) 


(2 - X)(l + k)C 

where A , B, and C are defined in equations (l8; and (20). 


(25) 


Equations (24) and (25) can he solved for 0^^ and C 3 /C L , after 
which C-^C^ and C 2 |C L can he found from equations (22) and (23). 

This work is best done with the numerical values for X, k, n, and 
Xq| c r for the particular wing under consideration. The procedure is 

illustrated by an example in a subsequent section. The relation 
between Xq/Cj. and x 1 / c corresponding to equation (13) for the 

pointed-tip wing is 


f o 2(1 + X + X 2 ) x_|_ 1 + X - 2X g 

c r = 3d + V ? 3(1 + \)(1 - *) 


(26) 


The spanwiBe load distribution corresponding to equation (15) for 
the pointed-tip wing can be obtained by substituting equations (22) 
and (23) into equations (21a) and (21b). The substitution into equa- 
tion (21b) produces only added complexity, and the result is not given 
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IT 


here. The substitution into equation (21a), however, gives the fol- 

X(l - k) 


lowing simplified equation for 0 § u ^ 1 - (k + n)(l - l) : 


cc l a 2(1 + X) 
Or Cq, « 


2(1 - X) (l - X 2 ) 

jt 


C4I 


r 2 - (1 - X)^t cj3 


(27) 


Unlike the spanwlse load distribution for the pointed-tip wing, 
the load distribution for the finite -taper wing cannot be conqpared .with 
an elliptical distribution for all wings but must be compared separately 
for each example investigated because of the form of equations ( 27 ) 
and (21b). The elliptical distribution is still given by equation (10). 


' NUMERICAL COMPUTATIONS 
Development of Form for Computation 


After numerical values have been found for the constants 0-^/0^, 

C 2/ C L, G 3/ C l' 81114 C^/Cl* these values are used with equations (4), 

or rather with the numbers in table I computed from these equations, 
to find the z-displacement corresponding to each component of the pres- 
sure distribution. The four displacements are then added to produce 
the fi n al shape of the warped wing. In principle this process is 
straightforward so that in practice it may be reduced to routine compu- 
tation. A form suitable for such computation is now developed. The 
particular form presented is one such that, at a given spanwlse 
station o, the z -ordinate as a fraction of the local chord c is given 
as a function of x'/c, the fractional distance behind the leading 
edge of the local chord. As a typical example, the z -ordinate corres- 
ponding to the second term of equation (2) is considered. From equa- 
tion (46) the following relation is written: 


m 



c 


Cg mc r 

R 2 n~ — T" 
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The following geometrical relations are easily verified: 


x 



, n x* -ffCi - x) 

1 - cr(l - X)j — + $ _ fe 7 - 


4 - 1 - ff(1 - x) 

mc r i . ]c 
s . " 1 - X 


so that 


Ha 


Cl c * c L 

Now If Aq Is defined as 


)[l - g(l - X)] 

LJ 

k' . cr(l - X) 1 

1 - X j 

|c + (1 - k) [l - o(l - xj] 


A 0(1 - X) 

*0 E (r- - x)j 


then 


_o_ (1 - fc)[l - g(i - X)~ 

Ao " 1 - X 

and. It can he shown that 

Ao x^ o(l - X) x' 

r c XT - k)Jl - o(l - XJJ " c + *o 

so that 


(28) 


m z 2 c 2 o /aA 2 

C L c " Rs C L Aq \r / 


(29) 


Thus, at a particular spanwiBe station o on a given wing, z 2 /c Is 

a function only of r (since R as a function of r is Ijnown from 
table I). From equation (28), x’/c is also a function of r; thus. 



( 30 ) 
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so that zg| c can be calculated for various values of x‘/c. By 

experience , it has been found that the most satisfactory way of choosing 
values of r to use in the computations is to plot, for the particular 
wing under consideration, x'/c against r (for the various values 
of cr) from equation (30). Only values of x'/c between 0 and 1 
Eire of interest. This fact determines, for each cr, the range of r 
to be used, and from the values of r used as arguments in table I 
those which give a satisfactory distribution of points along the chord 
are chosen. Relations similar to equation (29) for the other 
z -components can be derived; these relations are 



> (31) 


m z h _ CWff\ 2 AoY 

C T c C T lA j \r/ 

L L \ O/ N ' 

» — 

The amount of warp is seen to be directly proportional to C L and 
inversely proportional to m (for a given n) . 


y For some purposes, the wing shape can be more conveniently 
expressed in terms of z/c r and x/c r rather than in terms of z/c 

and x'/c. This conversion is easily made by use of the following 
equations: 


4 - [l - a(l - xj] f 


. (32) 


The foregoing equations are 
explanatory. This form has 
some of which are discussed 


embodied in table H, which is self- 
been used to compute several examples, 
in the next section. 
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Numerical Examples 


Example I. - The first example chosen has the following character- 
istics: 

n = 0.6 
k = 0.6 ; 

X = 0 , 

C* = 0 at 2L = 0.25 
c 

This set of characteristics represents a sweptback wing, tapered to a 
point at the tips, with the center of pressure a little more than 
50 percent of' the mean aerodynamic chord ahead of the location for the 
corresponding flat wing. The lift coefficient, sveepback angle of the 
leading edge, and Mach number are not specified} the final amount of 
warp is directly proportional to the lift coefficient, and any combi- 
nation of sweepback angle and Mach number that gives n = 0.6 may be 
chosen. — ■ ■ 


The four constants are found from equations (12) to have the 
following values: .... .. .14 ■ 

■ C-T 

= 4.4530 : 

C L ... - • 

7# * -19.0819 ‘ 

C L - 

C 3 

7T- - 15.9021 

C L • 
c 4 

^ = 0.9071 
°L 

Several values of a. are selected and, for these values, plots of 
x*/c against r are made from equation (30). Slide-mile accuracy 
is sufficient for these calculations, and only a few points need be 
taken to define the curves. The resulting curves for example I are 
shown in figure 7« These curves are used to' pick values of r from 
table I. The corresponding values of R-^, Rg, R^, and R^ from 

table I are entered in table II, together with the other necessary 
data, and the indicated computations are carried out. . 
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The results of the computations ere shown in figure 8, in which 
the ordinates are given as fractions of the local chord and the origin 
of the axes is at the trading edge of the local chord for each value 
of a. Several features of the wing are evident from this form of 
presentation: namely, the reflex curvature of the airfoil sections 

near the center of the wing (the angle of attack is infinite at the 
center line), the disappearance of this reflex curvature at outboard 
sections, the relative twist between inboard and outboard sections, 
and the (variable) dihedral. A better picture of the actual wing is 
obtained by plotting the results as in figure 9* In. order to give more 
physical maaning to 1 the picture, the results have been plotted for a 
lift coefficient of 0.2 and a leading -edge sweep of 60 ° (m = 0.577)* 

This last value thus corresponds to a Mach number of 1.44 since n 1 b 
equal to 0.6. 'There 'are two points worth mentioning with regard to 
figure 9. The first is that, within the accuracy of the linearized 
theory used in this paper, an arbitrary z(cr) may be added to the 
vertical ordinates without changing the aerodynamic characteristics of 
the wing. As pointed out in reference 1, this procedure is permissible 
so long as the resulting wing does not lie far from the z « 0 plane 
(that is, modification of the wing shape by addition of a set of ordi- 
nates which depends only on a (not on x) may be practiced in moder- 
ation) . The practical significance of this point is that the wing 
shape may be modified by this procedure to simplify the problem of 
locating spars. The other point is that for most configurations the 
inboard stations of the wing, which are those having the largest warp, 
are buried within a fuselage and, therefore, present no structural 
problems. (The effect of the fuselage on the aerodynamic character- 
istics is discussed subsequently. ) This particular examp le and pre- 
sumably others not too extreme should therefore be quite practical to 
build. 

The pressure distribution for this example is shown in figure 10. 
Because of the far-forward specified location of the center of pressure, 
part of the wing carries negative lift. The spanwise distribution of 
load is that shown already in figure 6. 

Example II .- The second illustrative example has the following 
characteristics : 


n = 0.8 
k = 0 
X = 0 

C m “ 0 at | r- 0.50 
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The values k = 0 and \ - 0 characterize a triangular wing. The 
center of pressure is at the same point as the center of pressure of 
the corresponding flat wing. The purpose of the present design is to 
show the kind of warp that might produce a wing with essentially the 
same center of pressure and spanwise load distribution as the flat 
triangle but without the steep pressure gradients that are known to 
promote leading-edge separation on the flat triangle, at leaBt at low 
Reynolds numbers. A constant-pressure triangular wing, of course, has 
the same center-of -pressure location as a flat triangle and has no 
adverse pressure gradients, but the spanwise load distribution of such 
a wing is triangular rather than elliptical. 

The method of computation is much the same as that used in the 
previous example. The principal difference is that equations (ll), 
rather than equations (12), may be used to find the following values 
for the four constants: 

Cl 

tF - 1.0907 


- -0.9082 

KC. 

c, 

=2 - 1.0907 

°L 

- 0.9071 
°L 


The results of computations made with tables I and II are presented in 
figure 11. These plots clearly depict a wing the main part of which 
is almost flat and which has a turned-down leading edge, a small twist 
from root to tip, and almost constant dihedral, angle alopg the span. 

For ease of manufacture, this dihedral can be removed without affecting 
the aerodynamic characteristic s very much. (See the discussion of this 
point under example I. ) The pressure distribution is shown in figure 12; 
the chordwlse gradients of pressure are not large. The spanwise load 
distribution is again given by figure 6. 
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Ryfijrrple III .- The third example chosen is a wing with finite taper, 
characterized by the following conditions: 

n «* 0.7 

k = 0.6 

X = OA 

C_ = 0 at = 0.25 

m c 

These values are substituted into equations (18), (20), and (26) to 
obtain the following values: 


— - 0.82857 
c r 

A = 0.05861 
B - 0.12308 

c - 0.20986 

Substitution of these values into equation (24) gives 

- -0.l8l4 

c L 

Substitution of this value into equation (25) gives 

r Cq 

p 3 - - 4.7438 
°L 

Th^ remaining constants are obtained by substitution of this value of 
into equations (22) and ( 23 )} thus. 

Si. . -5.2614 

C L 

^ - 2.6451 

t 

From this point the method of calculation is the same as that used in 
the two previous examples: suitable values of r are chosen, and the 

form of table II is followed to arrive at values for the wing ordinates. 
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The resulting wing shape is shown in figure. 13 together" with the pres- 
sure distribution; the results have been plotted for a lift coefficient - 
of 0.4 and a leading-edge sweep of. about 59° (m 0 0.6), corresponding 
to a Mach number of 1.54. The center-of -pressure location shown in 
the figure for the flat wing was found from reference 5* The wing 
shape is not extreme, and the previous remarks concerning the removal 
of the dihedral angle apply equally well to this case so that the wing 
can be built feasibly. The pressure plot shows the result of the 
assumption regarding the pressure in the tip region. 

# a 

The following spanwlse load distributions are found from equa- 


« 0.891 - 0.502 a 2 + 0a09ff3 


2.82 - 2.44c - 0.49c 2 + 0.11c 3 

These equations are plotted in figure l4, which also shows an ellipti- 
cal load distribution for comparison. The load distribution for the 
example being discussed is a fair approximation to the "ellipse so that 
no large drag Increase relative to the flat wing is to be expected as 
a result of the specified forward location of the center of pressure. 
As a matter of fact, the spanwise load distribution of the flat wing 
is itself not elliptical, so that the drag of the warped wing might 
well be less than that of the flat wing. ' 


NOTES ON PRACTICAL APPLICATION 


tions (27) and (21b): 
For 0 £ 0'S 0.795, 

cc z 

c r°L 

For 0.795 £ c £ 1, 
cc z 

c r C L = 


R an ge of applicability .- The method described in the preceding 
sections is directly applicable to wing plan forms of the types shown 
in figure 1. The locations of the various Mach lines shown in the fig- 
ure relative to the leading and trailing edges and relative to the 
center line are significant. The leading. edge must-be subsonic and the 
trailing edge supersonic; these conditions are expressed by the 
following inequality: 

l£n£|k| (33) 




TfiTr-,. 
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For the case of plan forms with finite tips (X ^ 0), the Mach line from 
the leading edge of one tip must not cross over to the opposite wing 
panel. This condition is expressed by the following inequality: 



In addition, because of the approximate nature of the assumption 
regarding the pressure in the tip region, cases in which the tip region 
covers a large part of the wing should be viewed with caution. 

Computing time .- The exact time required to compute a given example 
depends on such factors as whether X or k or both are equal to zero 
and the number of points taken to define the wing surface . The fol- 
lowing time estimates are given as representative of those required by 
using a manually operated calculating machine. To calculate the con- 
stants, l/2 to 1 hour is required, and to calculate eight spanwiBe 
stations, with 14 points along the chord at each station, 8 to 12 hours. 

Body effect .- In the derivation of the present method, the wing 
has been considered as Isolated) whereas in practice it is usually 
mounted on a body, on which may also be mounted a tail. The available 
information, both theoretical and experimental, is not yet sufficient 
to allow an accurate quantitative prediction of the effect of the body 
for the general case. (See references 6 and 7 for a discussion of the 
problem. ) Some qualitative estimates can be made, and by reference to 
whatever experimental data may be available for configurations resembling 
the particular example under consideration, rough quantitative correc- 
tions can be applied for the effect of the body. If the wing Is mounted 
on the body so that the chord line at the .Juncture Is parallel to the 
body center line, then the lift of the combination when the wing is at 
Its design position with respect to the free stream will probably be 
close to the sum of the lifts of the isolated wing and, the Isolated 
body. If, however, the two are connected so that when the wing Is at 
its design position with respect to the free stream the body Is at 
zero angle of attack, then the lift of the combination will probably 
be somewhat less than the design lift of the wing. In all cases, the 
center of pressure will probably be somewhat more rearward than that 
calculated when only the isolated components are considered. If these 
foregoing statements, which are obviously conjectural in nature, are 
accepted, then some allowance can be made for the effect of the body 
by adjusting the design conditions of the wing. The body- interference 
problem Is neither different nor more serious for the warped wings 
considered in this paper than for conventional flat wings, and all the 
preceding remarks apply equally well (or poorly, as the reader may 
judge for himself) to both types of wings. 
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Off -design operation . - In the course of a flight , the wing may he 
required to fly at the design Mach number at; attitudes other than that 
for -which it was designed. Within the limits of the linearized theory 
used in the analysis, the principle of superposition applies. The lift 
(and pitching moment) of the warped wing at an attitude different from 
the design condition is therefore simply the' design lift and pitching 
moment plus (or minus) the lift and pitching moment of a flat wing of 
the same plan form at an angle of attack equal to the angular deviation 
of the warped wing from its design attitude. | When the wing is required 
to operate at Mach numbers other than the design value, however, no 
simple method is available for estimating the change in aerodynamic 
characteristics, and even to calculate the properties by the use of the 
linearized theory is a practicably impossible Job. An experimental test 
Is the only way to find the answer. 


Applicability to other problems .- Although the derivation of the 
complete method has been limited to wings of" the types^shown In fig- 
ure 1, with approximately elliptical span loads, the basic results 
presented in equation (V) and table I are applicable to other wings as 
well. For example, a derivation similar to that presented in this 
paper could be made for sweptback wings with. cross-stream tips, such 
as that shown in figure 9 of reference 3* It is also conceivable that 
in some cases the shape of the spanwise -load-distribution curve mi edit 
be determined by some condition other than that of low Induced drag. 
The information presented in equation (4) and table I could be applied 
to such cases. 


As an example of an application of the basic data of equation (4) 
to a problem of a type different from that discussed in the section 
entitled "Numerical Examples", the design of a triangular wing with 
approximately elliptical loading in both the spanwise and the chord- 
wise directions is discussed. For convenience, this wing is called 
example IV. (in reference 8, Jones has shown that, for a lifting 
surface of narrow proportions lying near the center of the Mach cane, 
the mini mum value of the drag due to lift is achieved when both the 
spanwise and the chordwise loadings are elliptical.) 


The 
tlon (2) 


chordwise load distribution Is found from integration of equa- 
to be given by the following equation: 


Local lift ^1 x 
Total lift = Cl c r 


c 2 1 c 3\ /x \ 2 1 £4 

ft + 2 ft/^rj + 3 C L 


+ 


3 


(35) 
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If, as in the previous examples, the conditions of equation (ll) are 
applied to equation ( 9 ), then the spanwise load distribution is given 
"by equation ( 15 ), and the values of ^2j^L> and a;re 

those of equation (l4) : namely, 

£l _ £3 

C L ” °L 



Substitution of these values into equation (35) gives the following 
equation for the chordwise load distribution: 


Local lift £3 _x_ (k 3 £i\/V \ 2 . L 16VjlV 
Total lift C L c r + \jt " 2 CjJ^J “ 3nJ\c r J 


The 



chordwise load is now specified to be zero at the trailing edge 
= lV This procedure gives the following value for C 3 /c^: 



k 


_ 8 _ 

3* 


and the chordwise load becomes 


Local lift _ f. _8_\ _x_ 
Total lift _ \ “ 3*y c r 



(36) 


The load distribution given by equation ( 36 ) is compared with an 
ellipse in figure 15 • The spanwise load distribution is also repeated 
from figure 6 for the sake of easy comparison. 


The w ing shape is readily calculated from tables I and II and is 
shown in figure l6 for C L = 0.2, M = 1.2, and n =■ 0 . 3 . The drag 

coefficient for these conditions has been found by graphical integration 
to be approximately Cp = 0.0081. The drag coefficient of a flat wing 
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at the same conditions is 0.0091 if full leading -edge suction is 
assumed or 0 . 0l6l if no leading-edge suction: is assumed (no leading- 
edge suction has been assumed for the warped' wing) . 


CONCLUDING REMARKS : 


A method has been presented for designing a sweptback wing to 
have certain specified flight characteristics at supersonic speeds. 

For example, a wing of given plan form, operating at a given super- 
sonic Mach number, may be designed to have a specified lift coefficient, 
a specified center of pressure, and a nearly elliptical spanwise load 
distribution. As an aid in the calculations required for any specific 
case, certain basic data and a computational form are presented as 
tables. The procedure is illustrated by several examples. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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Figure 1-- Flan forms to Trial ch the method Is applicable. 
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Figure 2.- Axle system for pressure distribution, k 5 jjj-j 0 ts 
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Figure 4.- Concluded. 
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Figure 9-- Wing shape for exnmplft I with ordinates expressed as fractions 
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Figure 10.- Pressure distribution for example I. -> 0.2; nt ■ 0.577; 

M - 1.44. 
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Figure lft>- Spanvise load distribution for example III_ compared with 

elliptical load distribution. 
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(a) Chordwise load. 



(b) Spanvise load. 

Figure 15»- Load distribution for example XV. 
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